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ABSTRACT 


This thesis presents an operational Lanchester-type model 
of small-unit amphibious operations. This relatively simple 
model has been developed to demonstrate the basics of model 
building to the beginning student interested in amphibious 
warfare. The model is a time sequenced, deterministic, force- 
on-force combat model that is implemented on a digital computer. 
A brief discussion of considerations for modeling amphibious 
operations is given. The details of the model are presented 
for a specific amphibious-warfare scenario. Additionally, a 
computer terrain-contour-line plotting program is provided 
to assist the combat modeler to fit a parameterized-terrain 


to real terrain. 
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I. INTRODUCTION 


During the past several decades, combat models have been 
widely used to support military decisions. As the art of 
combat modeling becomes more advanced, combat modelers are 
continuously building more and more complicated models. To 
the beginning modeler, the ability to understand how those 
models operate is difficult, if not impossible. It is the 
purpose of this thesis to develop a simple amphibious-operations 
model that will demonstrate the basics of model building to the 
beginning student interested in amphibious warfare. In the 
broadest sense, an amphibious operation is a combined-arms 
operation which includes all forms of combat--land, air, sea. 
This thesis will limit itself to the small-unit amphibious 
operation. 

This project started with two basic models: one was de- 
veloped as the auxiliary model for the evaluation of design 
and employment alternatives for the LVA (Landing Vehicle 
Assault) in the thesis of David Larkin Chadwick (September 
1978). The other is the Smoler-Mills model which was developed 
in the thesis of Josef Smoler (September 1979) and enriched 
in the thesis of Glen Mills (September 1980). 

In Chapter II, general considerations for modeling amphi- 
bious operations are briefly discussed. Then, a small-scale 
computer-based Lanchester-type amphibious-operations model is 


presented, including analytical details of the algorithms used 








to represent each of the combat processes considered. Although 
the model was developed for a specific scenario, it is suffi- 
ciently general in design so that it can be adapted to other 
small-scale amphibious~operations scenarios with only rela- 


tively few modifications. 





II. GENERAL CONSIDERATIONS FOR MODELING 
AMPHIBIOUS OPERATIONS 
A. CHARACTERISTICS OF AMPHIBIOUS OPERATIONS 

In order to model an amphibious operation, it is first of 
all necessary to understand what is going on in a real amphibi- 
ous operation. Only after one knows the details of what is 
happening in such a complex combat operation, can one begin 
to sift out the cluttering details and make valid simplifying 
assumptions to come up with a tractable model. 

One of the key characteristics that serves to distinguish 
amphibious operations from other types of military operations 
is that a complete military force must be transferred ashore 
in an orderly manner under the constant pressure of actual or 
potential attack from hostile forces. Because the over-all 
amphibious assault requires precise and timely execution, the 
various component operations must be carried out in a planned 
sequence (especially early in the assault) according to a 
strict schedule. This sequence and these schedules, however, 
must be sufficiently flexible to permit rapid changes in line 
with unexpected development afloat and ashore. 

The notable success of amphibious operations during and 
after World War II is testimony to the fact that, with proper 
planning and organization, this dual problem can be solved. 
In today's environment, it is well known that the modern battle- 
field will be dominated by highly lethal weapons. This has 


raised serious questions about the survivability of amphibious 





forces. On the other hand, the long-range, high-speed assault 
potentially gives one the capability to launch assaults from 
far out to sea, land at times and places of one's choosing, 
and carry more firepower to accomplish the amphibious assault 
with greater safety for ships and men. Use of some type of 
combat model is the only way to explore such issues today. 

In order to build such a combat model of an amphibious opera- 
tion, it is necessary to develop and consider detailed and 
specific information on individual tactical and support ele- 
ments of the landing force, on the size, numbers, and charac- 
teristics of the equipments of these elements, and on the 
sequence of movement of these elements. 

The amphibious operation is a combined operation, the 
entire spectrum of activities involved in an amphibious 
operation includes: 

- pre-assault bombardment by ships and aircraft 
- sea mine clearance 
- attack on ships by enemy aircraft and cruise missile 
- ship-to-shore movement 
- surface assault landing 
- helicopter operation 
- ground combat between maneuver units 
- artillery and naval gunfire support 
- tactical aircraft support 
- mine warfare (sea and land). 
While an amphibious operation is one of the most complex of 


all military operations, defending against it is even more 
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complex: it is absolutely impossible for an enemy to defend 
all coastal areas at all times. The flexibility to conduct 
helicopter-borne vertical assault and surface-borne assault 
simultaneously will grealy enhance the complexity of defending 
against an amphibious assault. 

Airborne troops and supplies were valuable during the 
Second War, and further developments in that direction are 
under way. But whether or not modern airborne tactics and 
techniques have supplanted (in a practical sense) seaborne 
assaults (such as those used from 1942 through 1945 in the 
Pacific and elsewhere), it should be noted that the military 
problem of landing forces on shores held by an enemy remains. 
The emphasis in the future will most likely continue to be on 
having the ability to project forces from the sea onto a hos- 


tile shore and to hold such a beachhead. 


B. MODELING APPROACH 

All models of military operations must abstract from the 
real world. Since it is obvious that an engagement between 
modern military forces is a very complicated process, one has 
to abstract, aggregate, and interpolate in order to scale a 
combat process down to manageable size for modeling purposes. 
A variety of modeling approaches are available. These range 
from simple Lanchester-type models to highly complicated, com- 
puter assisted, high-resolution simulations in which the 
actions of each individual combatant are traced through a com- 


bat scenario second by second. Between these two extremes are 
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other approaches covering the whole spectrum of land combat, 
from one-on-one duels to theater-level models covering huge 
geographical dimensions. 

There are basically four different types of combat models: 
war games, analytical (or mathematical), simulations and some 
combination of these first three types. According to Bonder 
[Ref. 2], war games are not a feasible mechanism for analyzing 
a broad spectrum of system alternatives in a responsive manner 
to meet a planning cycle requirement. However, they are diag- 
nostic in the sense that they reveal problems that need to be 
resolved with future systems, and are a viable mechanism for 
training decision makers. Analytical models seek to describe 
the combat process mathematically. They simplify the conduct 
of sensitivity analysis and provide an increased ease in inter- 
preting results, since the dynamics of the combat process are 
contained in readily examined equations. Analytical models 
of any degree of complexity usually do not yield convenient 
analytical solutions but require numerical approximation 
methods. Simulation is the most widely used technique in 
Military system analysis. Simulation can generally produce 
very useful data, which are needed for further analysis, and 
sometimes for planning itself. However, the large amount of 
detail contained in most Monte Carlo simulations makes it diffi- 
cult to use as the sole vehicle to single out those systems 
capabilities, tactics, and environmental conditions which 
Significantly contribute to or delimit the system's effective- 


ness. Since, as we have seen above, no one type of combat model 
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is unconditionally preferred to another, it is proposed that 
a combat model should be selected or designed based on a 
specific scenario and upon analytical requirements. 

In most cases, detailed models are more credible to deci- 
Sion makers. However, for many people such detailed models 
of large-scale combat operations are far too complicated to 
be understood, require too much input data, and (in general) 
are not responsive enough. When one looks at computer storage 
and run time requirements for even the smallest high resolution 
model, it is easy to see why a high resolution model of a 
corps or theater is presently impractical, and is likely to 
remain so. In order to avoid the complexity of the large-scale 
model and to better understand land combat there is a growing 
trend among analysts to combine small unit and large unit 
models in such a way that the output data of a high-resolution 
small unit combat model is used as the input data for a low- 
resolution large unit model. The obvious drawback of this 
hierarchical-modeling approach is that any errors in the small 
unit models will be carried through, and possibly multiplied, 
as the process proceeds from model to model. In the large 
units the emphasis has been away from simulation and towards 
detailed Lanchester-based models. 

So far the emphasis has been adding more and more detail 
to the high resolution models so as to pick up as many inter- 
actions as possible. No matter how much detail is added to 
the small-unit simulation, it seems impossible that reality 


will ever be matched exactly. With this in mind, it is proposed 
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that a well-constructed Lanchester-based model of small unit 
engagements could give results that are just as valid as the 


results of a high-resolution simulation. 
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Til НЕ MODEL 


A. GENERAL 
1. The Scenario 

The scenario considers an amphibious-landing team, 
consisting of reconnaissance, a light infantry unit, and 
landing-assault vehicles. This team is part of an Amphibious 
Task Force (ATF), and it disembaxrks from ships that are on 
station over the horizon from the selected landing site. The 
assault vehicles, after transmitting from the amphibious 
shipping to the designated area for the landing formation, 
form into conventional landing waves at a distance offshore 
which is greater than the effective range of the direct-fire 
weapon systems of the shore-defense force. During the ship- 
to-shore movement the defender's anti-tank guided missile and 
improved gunnery system respond to the landing. Naval gun-fire 
ships provide fire support for the assault team during the 
ship-to-shore movement and the initial stages of landing. 

As the assault vehicles reach the beach, they (together 
with the assault vehicles and any weapons landed by landing 
vehicles) launch an attack on the enemy shore defense posi- 
tions. The defenders occupying those positions fight until 
their losses exceed a maximum permissible amount. The attack- 
ing force, however, continues the assault irregardless of 
losses incurred. Once the shore assault has been completed, 


the landing force with tactical mobility moves inland to 
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carry out the tasks, while the enemy prepares to mount a 
counter attack. 

The attacker may engage the advance force of the de- 
fender's initial counter attack force on the way to move in- 
land. The advantage will likely go to that force which has 
gained the initiative (i.e., the landing force) provided it 
can maintain its momentum. 

2. General Description of the Model 

The model developed in this thesis is a time sequenced, 
deterministic, force-on-force computerized model, coded in 
FORTRAN. The model conducts the battle in uniform time steps 
of 10 seconds each. Figure 1 shows the general scheme for the 
sequence of events in the model. The model simulates two main 
phases in the amphibious operations: (I) the amphibious-assault 
phase, and (II) the subsequent ground attack. The framework 
and the logical interrelationships of these two phases will 


be discussed in the following subsections. 


B. THE AMPHIBIOUS ASSAULT PHASE 
1. General 

In this phase the model considers attrition between 
the shore-defense force and the landing-assault force during 
its water-borne movement and subsequent assault to shore. 
The model aggregates the various actual combat organizations 
involved in the waterborne phase of the amphibious operation 
into several homogeneous combat units. Each of these units 


is characterized by certain relative offensive and defensive 
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READ INPUT DATA 
SHIP-TO-SHORE MOVEMENT 
SHORE ASSAULT 





REPEAT EACH 
TIME STEP 









UPDATE LOSSES 
PRINT OUTPUT 
ALLOCATE LANDING FORCE TO THE GROUND COMBAT 


MOVEMENT 


|GROUND COMBAT 
REPEAT EACH 
TIME STEP 

UPDATE LOSSES 


PRINT OUTPUT 


Figure 1. General Scheme for the Amphibious 
Operation Model 
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Capabilities. The following table illustrates the combat 
organizations that were explicitly modeled. The combat 


strength of each unit was represented by the state variables 


indicated. 
combat organization state variable 
Shore Defense--TANK assets DT 
Shore Defense--ATGM assets DS 
Incoming assault waves of LVA т = 12...05 


representing waves l through 5 

A cumulative combat force comprised TLF 
of those Marine ground units 

which have arrived at the beach and 

have debarked the LVA 


Fire Support Assets of The ATFFS 
Amphibious Task Force 


The initial strength in each of the above units is input 
data to the model. 

The schematic of the method of employment for the LVA 
in the ship-to-shore phase of an amphibious assault is shown 
in Figure 2. It is assumed that the conventional landing 
formation composed of waves of landing vehicles will be used 
as prescribed by current doctrine. The movement of assault 
vehicles to the beach is simulated using a time step approach. 
At each time increment the positions of vehicles are updated. 

The tactical interrelationships which exist between 
various combat units are illustrated in Figure 3. Assuming 


that in such a future amphibious operation the attrition of 
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Figure 2. Concept of Ship to Shore Movement 
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incoming landing vehicles would be dominated by the effects 
of shore defense direct-fire weapon systems (specifically, 
tank and anti-tank guided missile (ATGM) assets), the model 
essentially omits the effect of the defensive indirect fire 
capabilities. 

2. Attrition Process 

The model represents the attrition of all combatant 
units as a deterministic process. The primary consideration 
in the ship-to-shore movement of incoming waves of assault 
vehicle is the attrition effects upon those waves due to the 
direct fire weapon assets ashore. The attrition of each wave 
utilizes Lanchester "aimed-fire" equations with variable attri- 
tion-rate coefficients. 

The classical Lanchester hypothesis for aimed-fire 
attrition (combat under "modern condition") is that the casualty 
rate of a unit is proportional to the "size" of the opposing 
forces. If the unit "A" is being engaged by "D", this may be 
expressed by the differential equation: 

= = - BETA), X D 
The proportionality constant BETA), is called the Lanchester 
attrition-rate coefficient. It is assumed that this func- 
tional relationship holds for each (firing unit, target unit) 
pairing over a small time interval dt. The problem then is to 
determine numerical values for the Lanchester attrition-rate 


coefficients. In this model, these coefficients were expressed 
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as the product of the rate of fire (ROF) and the kill proba- 
bility per round (Р). Thus, 


BETA DA = ROF pA x PX) ΠΑ 


The rate of fire (ROF) can be expressed as the reciprocal of 
TBF (Time Between Firings) which can be evaluated by 


TARGET RANGE 


ΠΤ ALM-RELORD TÍME + TARGET SPEED + PROJECTIVE VELOCITY 


In determining the probability of a vehicle "KILL" per 
round, it is assumed that a hit by a large caliber projectile 
would constitute a "KILL" and the two defensive weapon sys- 
tems addressed would exhibit normal, uncorrelated horizontal 
and vertical errors. Then the single shot kill probability 


is given by 


a-u/co, 2 


1 exp(- =—) ах] 


/2π (-a-u)/o, 





P(k) = Г( 


b- 
( τυ. 


- ехр(- ἅτ) άγ] 


πὶ || 
"PE (-b-v) /o., 





where: 


a = semilength of a target 
b = semiwidth of a target 
u = horizontal aiming error 
у = vertical aiming error 
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round-to-round standard deviation in vertical 


Q 
il 


σ round-tround standard deviation in horizontal 


y 


Model functions RNG, HT and SPD are called upon within 
the model logic to generate the range, height and speed re- 
spectively for each assault wave as time is incremented through- 
out the course of the amphibious assault phase. This informa- 
tion and typical dispersion data (both mean and standard devia- 
tion for the tank and ATGM weapons) are then incorporated into 
the rate of fire and hit probability calculations. 

The Amphibious Task Force's fire support assets contri- 
bute significantly to the combat effectiveness of the shore 
defense units. Since it is assumed that the exact positions 
of the defensive units DT and DS emplaced on shore are unknown 
to the Amphibious Task Force and consequently the ATF fires 
into the general areas thought to contain the defensive units. 
The following Lanchester-type area-fire equations are applied 


to compute the attrition of DT and DS. 


арт _ 

<: = -(ALPHA,, x ATFFS) x DT 
dDS _ (ALPHA... x ATFFS) x DS 
dt DS 


The combat effectiveness of the ATF fire support assets is to 
be considered relatively constant during this segment of com- 
bat time. Thus the terms in parentheses on the right hand 


side of these equations are to be considered an input parameter. 
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Once a particular defensive unit has initiated its 
engagement of incoming waves it is considered that their 
fire "gives away" their exact locations. At this point it is 
assumed that the ATF fires will engage that defensive unit 
through the use of aimed-fire and the loss rate will be in 
accordance with the Lanchester hypothesis for aimed fire. That 


is, 


ат _ _ 
арт _ 
Pru = ВЕТА x ATFFS 


Again, the parameters on the right-hand sides of both these 
equations are provided as input. 

The casualty rates applied against the DT and DS by 
the Total Landed Force (TLF) are determined by means of the 
Lanchester aimed-fire attrition rate coefficients by the 


equations 


dt ^ "WBETAQ p-pT * TF pr 
арз _ _ 
=> ИВЕТА, „рс х ТЕС 


The computation of these WBETA coefficients is not 
performed within the model utilizing the detailed rate of fire 
and P(k) arguments described previously but is required as 


input. Although the defensive losses are considered significant, 
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a high level of complexity for computing these coefficients 
has not been incorporated into the model at this time. 

Figure 4 describes the schematic of the attrition 
process of the amphibious assault phase in the model. The 
attrition during each time step was computed using the Euler 
integration method to approximate Lanchester's force-on-force 
attrition differential equations. 

3. Fire Allocation 

Each weapon category was asSigned an engagement window 
as illustrated in Figure 3. Only those LVA located within 
these range windows could be fired upon by the shore defenders. 
A defensive weapon only engages the two closest incoming waves 
if more than two waves of LVA are at any time located within 
the weapon's engagement window. If only one wave of LVA is 
present in a weapon's engagement window, defensive fires of 
that particular weapon type will be distributed uniformly 
against the surviving LVA in that wave. 

If two waves of LVA are both contained within the 
engagement window, defensive fires of that particular weapon 
type will be distributed according to a tactical allocation 
submodel. A weighting factor (DEFWT) is utilized in estab- 
lishing the proportion of the total weapon strength to be 
allocated against the surviving LVAs in each of the two waves. 
As an example, if DEFWT(l) = 2 and DEFWT(2) = 1, then each 
surviving LVA in the closer of the two incoming waves would 


be allocated twice as much fire as surviving LVA in the seaward 
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———ы = 


wave. For the purpose of these examples, if waves 3 and 4 were 
both located within the tank engagement window, then the pro- 
portion of DT's fire allocated to surviving LVA in wave 3 


would be 


DEFWT (1) x WV(3) 


DEFWT(1) x WV(3) « DEFWT(2) x wv(4j * PT 


where WV(3) is the state variable for the current number of 
survivors in wave 3. 

As each assault wave arrives at the beach, the total 
surviving strength of that wave is transferred to the variable 
TLF (Total Landed Force). TLF represents a ground combat 
force equal to that transported by the number of LVA survivors 
having arrived ashore. Once established, TLF engages the two 
shore defensive units allocating its fires between the two 
defensive weapon categories in the same proportion as the 


number of surviving tanks and ATGM's, that is 


Тс = pr F pg * TL 
E DT 
pp s pr os S Ish 


4. Suppression 


The suppression effects of incoming fire upon each of 
the defensive units was considered a significant factor with 
respect to its effect on the survivability of the incoming 


assault waves of LVA. Generally, the effect of suppression 
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fire will hinder an individual from observing and firing at 
the enemy. 

It was assumed that suppression would degrade unit 
effectiveness by increasing the aim-reload time (ARTM) and 
round-to-round error standard deviation for each weapon 
system. Hypothesizing that attrition rate is the dominating 
variable, and therefore, a good indicator of the suppression 
level, ARTM and such round-to-round errors were assumed to be 
functions of the force's attrition rate. This is an area, 


however, requiring further study. Analytically, 


ARTM = ARTM (1 + GAMMA x DA) 
sup nonsup 
ERROR, D = ERROR, onsup ‘+ + DELTA x DA) 
where: 
DA = attrition rate for defensive unit due 


to the effect of AFTTS and TLF 


GAMMA = parameter representing relationship 
between DA and ARTM 


DELTA = parameter representing relationship 
between DA and error standard deviation 
This increase of ARTM and round-to-round error (expressed as 
a standard deviation) decreases the kill probability (Ρις) tor 
both defensive weapon categories. Parameter estimation would 
appear to be the largest problem. But, since determining 
these parameters in the model is beyond the scope of this 


thesis, these parameters GAMMA, DELTA are provided as input. 
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5. The termination of the Assault Phase 
It is assumed that if during the course of the amphibi- 
ous operation the shore defense forces suffer a cumulative 
loss in excess of 70% of their initial force strength, the 
remaining shore defense will try to withdraw, resulting in 


termination of the engagement. 


C. THE INITIATION OF GROUND ATTACK 

In the amphibious operation, the landing-force must sieze 
critically-important inland objectives as rapidly as possible 
before the defenders start to react to the landing. The deci- 
sion for the initiation of ground attack should be based upon 
the enemy threat and desired landing-force build-up ashore, 
among other factors. To model this decision rule, it is 
assumed that once the landing has begun, the landing-force 
commander will base his decision about initiation of ground 
attack primarily on the strength of the landing-froce ashore 
and the shore defender's strength. The criteria for the deci- 
sion should meet these two conditions: 

(1) The survived landing-force strength has to be 
greater than the minimum force required to carry 
out the ground attack. 

(2) The defender's strength must fall below the minimum 
required to continue coordinated shore defense before 
breakoff and retreating. 

These conditions are then checked after each time step. If all 


waves landed without reaching the above second conditions, it 
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is assumed that the next wave group will engage any leftover 
defenders. Thus, the decision to implement the ground attack 


is based on the size of the total landed force. 


D. THE GROUND ATTACK PHASE 

The attacking force which is composed of three subunits 
of three LVAs armed with TOW antitank missile system attacks 
along predetermined routes. The defending force is comprised 
of three subunits of three tanks in a static defense. 

The battle takes place on parameterized terrain which will 
be discussed later. The ground-attack process contains five 
main subprocesses: (1) movement, (2) detection, (3) fire- 
allocation, (4) attrition and (5) battle termination. The 
general flow of the ground attack phase is shown in Figure 5. 

l. Movement Process 

Every attack unit is advanced to the next interval 
along a predetermined route unless this unit is destroyed 
already or is in firing status. To use his own determined 
routes, the user is required to input the original location 
of each attacking subunit and from one to ten nodes he wishes 
each attacking subunit to move through. This information, 
along with vehicle speed, is used to calculate route intervals 
that move the attacking unit through each of the designated 
nodes. The straight line ground distance between the first 


two adjacent nodes, DIST, is calculated as 


DIST = ХЕ У 
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initiate ground attack 


update locations of attacking units 


check line-of-sight between any two opponent units 






fire conditicns allocate 
are satisfied targets 
update the accumulated detection probabilities 
repeat each 


- campute attrition 
print results for the current time step 


Figure 5. Flow Chart for the Ground Attack Phase 
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where: 
X = distance between two nodes in straight west- 
east direction 
Y = distance between two nodes in straight south 
to north direction 

The angle between the desired direction of movement 
and straight west to east movement, a, is then calculated. 
Utilizing these quantities, the distance desired to move 
during each time step (DST) and the distance to be moved in 
the X and Y direction (XLN,YLN) is computed. These distances 
(XLN,YLN) are then added to the coordinates of the previous 
interval endpoint to determine the coordinates of the next 
interval endpoint. This process is repeated between the next 
two nodes until the unit has traversed the entire route. 

2. Detection and Fire Allocation 

The detection phenomena are modeled in two ways: 
(I) non-firing detection, and (II) firing-detection. A non- 
firing detection can occur as a result of an observer's random 
search within his designated sector of responsibility. The 
probability (Ру) that an observer A looking at the direction 
which enables him to detect a target is computed by integrating 
the Limicon Function over limits are +15° from the primary 
direction the observer is looking. The Limicon Function, f(9), 


is the following probability density function: 


f(0) = А + B x cos(9) EU uo 0 
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where: 


a AAA —< 


assigned sector width/2 


- B x cos(D) 


>(sin (D) - D x cos (D)) 


primary direction observer is looking 


Assuming 30? field of view for any observer A target 


B might be seen only if the observer A is looking at the 


direction such that ANGRT « 0 « ANGLFT where: 


ANGLE 


ANGLFT 


ANGRT 


Thus 


the absolute value of the angle between 
the primary direction (IPRDIR) and the 
observer-target direction (OTANG) 
апае lo τε angle + 15 D 

D 1f angle + 15 > D 


angle - 15° 


pr (ANGRT « a < ANGLFT) 


ANGLET 
f £(a)da 


ANGRT 


Given that observer A is looking at the direction, the 


curve [Ref. 11] 


conditional detection rate (Хав) is determined by the regression 


The probability that unit j is detected by 


unit i at time Cat [P,, (t + At)] is computed according to: 


t+At 
- Í à (t)dt 
E 
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gs can be interpreted as the average fraction of unit i 
that detects unit j. 

The second method of detection played in the model is 
a so-called firing detection. This phenomena occurs when the 
following happens: if a firing location is within 1155 ος 
an observer's primary direction of observation when he is 
firing, he is assumed to be detected and is added to the 
observer's target list. In summary, the following conditions 
are necessary for unit j to be a target of unit i: 

(a) Line-of-sight must exist between unit i and unit j. 

(b) The range between the two units should be between 
maximum range and minimum range of the firer's weapon 
system. 

KG) Er - At) » O. 

The fire-allocation routine determines what fraction 
of each unit is allocated to fire each target in target list 
since it is assumed that each fire unit is not restricted to 
fire at one target. This fraction is determined as a function 
of the predetermined fire policy and Den The fire policy 


is as follows: 


# of 3 ОЕ unit i allocated to each target 


EJ st nd rd 
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The priority of a target is taken to be a function of range 
only. The fire allocation rule which is used in the model is 
documented in detail in Smoller's thesis [Ref. 9]. 

Dee Aer Lt Lon 

The attrition process in the ground attack phase uti- 

lizes Lanchester "aimed-fire" equation used with variable 
attrition coefficients. The calculation of the attrition 
coefficients is accomplished through the use of one of two 
optional methods. The first option uses the following Bonder- 
Farrell formula to compute the reciprocal of the expected 
time to kill. The coefficients, Αλαν the rate at which one 
firer of unit 1 kills unit j targets are computed according 
ЕО: 


Aij = LE ar 


where D is the expected time for one firer of unit i 


to kill one target of unit j. The ΕΝ is computed using: 


E (T; .) = t, + ti - t 


j + (ty t te) /P(K/H) 


h 


* (tQ * t) /P(h/m) x ((1-P(h/h)) /P (K/H) 


РВ) 


where: 


(t 
I 


time to acquire a target 


ct 
il 


time to fire first round following acquisition 
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Ë = time to fire following a hit 


h 

ЕЗ ТЕ time to fire following a miss 

te = time of flight of a round 

P = probability of a first round hit 

P(h/h) = probability of a hit following a hit 
P(h/m) = probability of a hit following a miss 
P(k/h) = probability of a kill given a hit 


This formula holds for the conditions that the hit probability 
of any round depends only on the result of the previous round 
and no accumulated damage is considered. It is assumed that 
P(K/H) = 1.0 and P(h/m) = p(h/h) = P, thus reducing the equa- 


Eion to: 


+ Е. + (t, τε) (1-P) /P 


a l 


EE = t +t £ 


The second method, called the stochastic method, 
interprets the attrition rate coefficient, А; as a measure 
of the fighting ability of a unit which has a random phenomena 
affected by many different factors. It is assumed that the 
random fighting ability should be distributed between .3 and 
.8 with the majority of the unit being rated between .5 and 
.6. A "fitted" distribution to these assumed fighting levels 


which is devised by Mills is: 


AD. = -202 + 2U + „3; U ls a random Uniform (0,1) number 


1j 


The A,,'s are a realization of the random variable denoting a 


unit's initial fighting capability prior to the battle. Then, 
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during each time step, a new attrition rate coefficient for 


each unit is computed using the equation: 


0 
A, (1 г/г.) Or OR - r. 
А.. = 
1) 
0 FOr r «тү 
e — 
where: 

а = maximum effective range of a firer's weapon 
r = current range between firer and target 


Utilizing one of the above formulas to calculate А;.'5, the 
attrition during each time step was computed using the Euler- 
Cauchy differencing equations to approximate Lanchester's 
force-on-force attrition differential equations. 
4. Termination of Ground Attack 
The ground attack is terminated when either: 
(1) One of the two opponent forces is annihilated; 
(2) A distance between each attacking subunit and each 
defensive subunit which is still engaging becomes 
"too close"; 
(3) Any attacking subunit passes by the flanks of the 
forward most defensive subunit still in the battle. 
The criteria for being "too close" is user input. This allows 
for flexibility of breakpoint distance for various weapon 


systems on. the battlefield. 
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E. THE PARAMETRIC TERRAIN 

The terrain affects a great deal on detection, mobility, 
tactics, and intervisibility between weapon systems in ground 
combat environments. In the model, the battle is simulated 
on 3 x 4 Km piece of terrain represented as a part of the 
coastal area east of the Korean Peninsula. It is important 
to have a terrain representation to emulate actual terrain 
areas. The model uses the parametric terrain representation 
method which was proposed by Chris Needle [Ref. 8]. The idea 
of the parameterized terrain is that the elevation of any hill 
mass can be represented by a bivariate normal density function. 
Mathematically, if f (X,Y) is a function giving the elevation 
of the I's hill masses at any X,Y map coordinates on the battle- 
field, the overall terrain elevation at X,Y is obtained as 


the positive maximum over all the hill masses, 


Z = £(X,Y) = maximum fr (X, Y) 
Е NHILLS 


where NHILLs is the total number of hill masses on the battle- 
field. Then, elevation data is used to compute the existence 
of line of sight between opposing forces which is a key element 
in detection process. The model uses the line-of-sight routine 
which was written by Prof. James K. Hartmann [Ref. 5]. 

In order to represent a piece of real terrain with parametric 
terrain, it is necessary to fit hill mass functions f (X,Y) to 


a contour map of the terrain to be modeled. The fitting 
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process can be done by comparing a computer generated contour 
map by varying the bivariate normal parameters to the original 
terrain map. The computer generated terrain map of the battle 
area is inclosed as Figure 6. Appendix D presents the program 
listing for plotting a contour map from hill mass functions. 
This program can be used for the user to fit a specific real 


terrain which he has in mind into the parameterized terrain. 
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IV. FINAL REMARKS 


In order to illustrate basic modeling points for under- 
standing and developing complex-operational amphibious-warfare 
models, a simplified model (specially tailored to a small-unit- 
amphibious-operation scenario) has been developed. Although 
several user options and varieties of modeling techniques 
have been incorporated into the model, the model should not 
be viewed as the final product. 

There are several areas of model enhancement and enrichment 
that should be considered in the future. Only the aggregated 
amphibious task force fire support is played explicitly in the 
amphibious assault phase. The effect of artillery, naval gun- 
fire and close air support during the various phases of amphibi- 
ous operations should be added to the model. Ammunition con- 
sumption and resupply are considered vital to the success of 
any military operation, and such aspects should be also added 
to the model. Enemy forces should be played in greater detail. 
The movement of reacting enemy force and their dynamic-defen- 
Sive-position selection are being considered for future model 
inclusion. A module which considers enemy reaction time, the 
terrain, and the tactical situation, and which dynamically 
determines which side is going to take the defensive role, as 
well as its defensive position, has been proposed to enhance 
the existing model. Inclusion of these features would create 
a more complicated model, adding realism but detracting from 


the current simple and transparent form. 
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The model currently simulates the ground combat on a 3 x 4 
Km piece of coastal terrain representing an area of Ham-Hung, 
Korea. In order to simulate combat on a user's specific ter- 
rain, he needs to fit a parameterized terrain model to the 
actual terrain of interest. Doing so using the terrain fitting 
method described in the preceding chapter is an extremely 
tediuous and time consuming process. The development of a 
more efficient terrain fitting technique will greatly enhance 
the model flexibility, and increase its utility in responding 
to interrogations from the real world. 

The interested reader may obtain the model program deck 
and the sample input data deck from Prof. James G. Taylor, 
Naval Postgraduate School. Since the data used for the model 
run is hypothetical and greatly simplified, the reader is 
cautioned not to draw any analytical conclusions from the 


putput. 
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APPENDIX A 


SAMPLE EXPECTED OUTPUT 


Amphibious Assault Summary 

AMPHIBIOUS ASSAULT INFORMATION 
INITIAL FORCE STRENGTH 

WAVE l 2 3 4 5 
LVA νυ. ο το 16ο 3.0 
DT = 3.0 DS = 1.0 

LVA ENGR SPECS 

SPDMAX SPDMIN HTMAX  HTMIN  WID 
10.30 3.50 1.770 0.60 3.53 


DEFENSIVE TACTICAL PARAMETERS 


RANGE AIM-RELOAD PROJECTILE 
MAX MIN TIME VELOCITY 
TANK 1500.0 15.00 350.00 
ATGM 2000.0 200.0 30.00 350.00 


DEFENSIVE TACTICAL ALLOCATION WEIGHTS: 
WAVE 1 = 2.00 WAVE 2 = 1.00 


DEFENSIVE FORCE ATTRITION COEFFICIENTS 


ALPHA*A BETA*A 
DT 0.00006 0.00070 
DS 0.00008 0.00090 
WBETA (1) = 0.00050 ‘WBETA(2) = 0.00070 


BREAKPOINT ASSUMPTION: 0.3* (TOTAL DEP FORCE) 
DEFENSIVE FORCE LEVEL FOR GROUND ATK 0.32 


DISPERSION DATA 


RANGE TSIGV RANGE TSIGH RANGE TMEANH 
25.0 0.0 25.0 0.0 25.0 0.0 
500.0 2.0 500.0 2.0 500.0 1:0 
1000.0 2.0 1000.0 3.0 1000.0 5.0 
2000.0 20.0 2000.0 20.0 2000.0 10.0 
5000.0 250 5000.0 25.0 5000.0 15.0 
10000.0 25.0 10000.0 25.0 10000.0 15.0 
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LI. 


EII. 


RANGE SSIGV RANGE SSIGH 
20 0.0 250 0.0 
250.0 5.0 250.0 5.0 
500.0 T3 50070 um 
1000.0 14.0 1000.0 14.0 
2500.0 15.5 2500.0 5-5 
5000.0 17.0 5000.0 17.0 
10000.0 20.0 10000.0 20.0 
Assault Phase Time Step Summary 
TIME = 815.0 SECONDS 
WAVE FORCE LEVEL STATUS LOST-PCT TSURV 
T ο το. 1 0 927 
2 1.7451 1 0.841 
3 4.5060 1 0.590 
4 9 5051 3 0.019 
5 3.0000 2 0.0 
TANK 0.0 1.000 19,79 
ATGM 0.0 1.000 0.0 
FINAL LVA SURVIVORS ASHORE = το 788 


GROUND ATTACK STARTS AFTER DEFENDER BROKE CONTACT 


GROUND ATK TIME = 825.0 


Initial Ground Combat Summary 


INITIAL GROUND COMBAT INFORMATION 


UNIT X Y FORCE LEVEL 
1 2000.0 1900.0 3.0 
2 1900.0 2400.0 3.0 
3 1500.0 2100.0 3.0 
4 3800.0 2700.0 3.0 
5 3800.0 2300.0 3.0 
6 3600.0 1700.0 320 


ATTRITION IS DETERMINISTIC 
ROUTES DETERMINED BY USER 
ATTACK VEHICLE SPEED IS 12.0 


BREAKPOINT DISTANCE IS 500.0 


DEFENDER WILL MOVE TO ALTERNATE POSITIONS 
ALTERNATE POSITIONS ARE: 


UNIT X M 
à 4500.0 3800.0 
2 4500.0 2700.0 
6 4600.0 1800.0 





Να“ 


TIME 


ATK KILL PROBABILITIES 


RANGE 

500 
1000 
1500 
2000 
2500 
3000 


DEF. 


RANGE 

500 
1000 
1500 
2000 
2500 
3000 


Ground Combat Time Step Summary 


UNIT 


Jt 
l 
2 
3 
4 
5 
6 
D 


E 
0.85 
0.80 
0% 75 
0.60 
0.45 
0.20 


KILL 
E 
0.60 
0.85 
0.80 
075 
0.60 
0.40 


1395 SECONDS 


X 
2420.8 
3664.1 
4397.4 
4500.0 
4500.0 
4600.0 


Y 
1984.2 
225320 
1742.2 
3800.0 
2700.0 
1800.0 


PHH PHM 
0.85 0.75 
0.80 0.75 
0.75 0.70 
0.65 0.60 
0.50 0.50 
0.20 020 

PROBABILITIES 

PHH PHM 
0.70 0.65 
0.90 0.85 
0.85 0.85 
0.80 0.75 
0.70 0.65 
0.45 0.40 


ISTANCE BETWEEN FORCES I 
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PKH 
0.70 
0.70 


0259 
0.35 
0.20 


PKH 
0-85 
0.90 
0.80 
0.70 
0.65 
0.50 


FORCE LEVEL 
0 


OO CLOSE. 


Ei toO to Oo rto 


LOST=PCT TARGETS 


1.000 
1.000 
9 059 
1.000 
0.055 
1.000 


ND OF BATTLE 





| 
1 


Amphibious Assault Input 


L 
0.30 


10. 


0 
0 


EI. 


3 


0.00006 


.0007 
.0005 
32 


50. 


Terrain 


46 

0. 
2000. 
1800. 
2000. 
2400. 
2450. 
2700. 
3200. 
4300. 
3750. 
4150. 
3200. 
4600. 
4800. 
2200. 
2400. 
3100. 


Data 


1100. 
2200. 
1900. 
1400. 
1700. 


1800 


1650. 
1300. 
1750. 
1600. 
2150. 
1700. 
1500. 
2600. 
2850. 
2700. 


APPENDIX B 


170. 
150. 
150. 
150. 
T30; 


140. 
130. 
150. 
150. 
130. 


170. 
170. 
150. 
T50: 


LISTING OF SAMPLE INPUTS 


e 


99979 
350. 
300. 


500. 
500. 


400. 
660. 
550. 
500. 
300. 
300. 
350. 
300. 
350. 


10000. 
10000. 
10000. 
5000. 

17% 


5000. 
by 


N) P= = IÑ Ко к CO CO CO C) NN) S М М Ñ OO 


Oy UN м м (л 


о ол лл 


10000. 


10000. 
20. 





2500. 2400. W50 


1 28943 

03 03 

0000.0 2500.0 0500. 4000. 
550. 3.0 3.0 

E 2 

2000.0 1900.0 

1900.0 2400.0 

1500.0 2100.0 


5000.0 2500.0 
4900.0 2150.0 
2200.0 1700.0 


4800.0 1750.0 
3800.0 2700.0 3.0 190 120 


2650. 2850. 150. 
2700. 2600. 150. 
3800. 2200. 150 
4500. 2600. 150 
3600. 2800. 150 
2700. 3300. 190 
3000. 3300. 170 
8150. 3750, 130. 
B750. 3200. 150% 
3800. 3800. 150. 
3600. 3600. 1550. 
4150. 3950. 170. 
1650. 2100. 150. 
2250. 2100. 180. 
4000. 2200. 150 
3900. 2200. 150 
0 0 0 0 0 
0 33 39 93 62 
0 0 0 0 0 
0 6 14 9 12 
101 
1 2 3 30 4 
6 32 33 7 TL 
8 9 10 11 33 
8 42 2 14 30 
16 17 18 19 20 
2 31 11 16 20 
46 20 21 22 12 
42 45 46 14 23 
26 14 25 26 27 
35 44 26 27 29 
40 
III. Ground Combat Input 


0 
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250% 


300. 


O UU O YN y UU O O O 


ω N FP N N ϱ 0 UN ο ο N F2 = = UC FP 
NONO 


сл 





ооо ооо 


ооо о с 


2300.0 3.0 190 120 
1700.0 3.0 180 120 


0 4 

3800.0 
2700.0 
1800.0 
95907590. 70 
1055 
το 05090 65 
65 0.60 0.55 
25090: 5090255 
2070720 0.20 
700,65 0.85 
-90 0.84 0.90 
8570783 0.50 
50: 0.750.710 
7/050, Go 0.605 
45 0.40 0.50 
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APPENDIX C 


DEFINITION OF VARIABLES IN COMPUTER PROGRAM 


1. The Amphibious Assault Phase 


CDSURV (I) 


CSURV (I) 


DA (I) 


051 


DS2 


DTl 


DT2 


DTIPH 


DTlROF 
DINIT 
IL(I) 
IWPN 


IWSTAT(I) 


Current strength of defensive force I 


l TANK 
2 ATGM 


I 
I 
Current strength of assault wave I 


Attrition rate for def. unit I due to the effects 
of ATFFS/TLF 


That portion of the DS unit assigned to engaging 
the closer of two multiple waves in the ATGM 
engagement window 


That portion of the Ds unit assigned to engaging 
the farther of two multiple waves in the ATGM 
engagement window 


That portion of the Dt unit assigned to engaging 
the closer of two multiple waves in the TANK 
engagement window 


That portion of the DT unit assigned to engaging 
the farther of two multiple waves in the TANK 
engagement window 


Hit probability of rounds fired by DTl against 
wave TENG(1) 


Rate of fire utilized by DTl against wave ТЕМС (1) 
Initial strength of def..force I 

When equal to l indicates the wave landed shore 
Weapon code: TANK = 1, ATGM = 2 

Current status of wave I 

- not engaging 

- landed 

- under fire by ATGM 


- under fire by TANK 
- under fire both ATGM and TANK 


m WN PO 
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GALF 


GATK 


GATM 


RD 


RKSURV (1) 
SA (1) 


SENG (1) 


SRNG (I) 
SSIGH 
SSIGV 


SWTS (I) 


TA (I) 


TBW 


TMEANH 
TMEANV 


TENG (I) 


TRNG(I) 
TSIGH 
TSIGV 
TSURV 


TWTS (I) 


WVINT(I) 


Denote whether the LF build-up is sufficient 
for the ground attack 
0 - not sufficient 
1 - sufficient 
Denote whether the LF initiated the ground attack 
0 - not started yet 
1 - started already 
Time when the ground attack started 


Distance offshore at which waves initiate 
transition 


Concatenation of CSURV and CDSURV 
Attrition rate for wave I due to ATGM 


The wave number of the closer of two waves in 
the ATGM engagement window 


Firing range to wave SENG(I) 
The std dev error in the horizontal for ATGM 
The std dev error in the vertical for ATGM 


The proportion of the total DS strength to be 
allowed to engaging SENG(I) 


Attrition rate for wave I due to TANK 


The interarrival time between waves arriving at 
the beach 


The bias error in the horizontal for TANK 
The bias error in the vertical for TANK 


The wave number of the closer of two waves in 
the tank engagement window 


the firing range to wave TENG(I) 

The std dev error in the horizontal for TANK 

The std dev error in the vertical for TANK 

Total number of surviving LVA at the current time 


The proportion of the total DT strength to be 
allowed to engaging TENG(I) 


Initial strength of wave I 


Sr 





2. The Ground Attack Phase 


ALPHA (1) 
APOA(I,J) 


AVSP 


BREAK 
DISMAX 
DIST 
DST 


FL(I) 
FO(I) 
FO(I) 


IALT 


IC 


IDIR(I,J) 
ЕЕ (І) 


IMOVE 


IPRDIR (I) 


IRTE 


IS 


ISECWD(I) 


Initial attrition-rate coefficient for 
stochastic attrition module 


The average proportion of the ¿En attacker of 
unit i allocated to fire on unit i 


Average speed of moving attacking units 


Breakpoint distance between attacking units 
and defenders 


Maximum distance allowed between attacking units 
before the leading unit is delayed 


The straight-line distance between two movement 
nodes inputed by the user 


The distance in meters to be moved each time 
step by attacking units 


Force level of unit 1 
Force level of unit i 
Initial force level of unit i 
Denotes whether the user desires alternate 
defensive positions or not 
0 - yes 


1 - по 


Counts number of time units a defender has 
been moving 


Direction of zen interval in ne route 
Interval index for unit i 


Number of time units a defender is allowed for 
moving to an alternate position 


Primary direction of movement for unit i 
Denotes whether user wants to input routes or 
NOE 

0 - program determined routes 

1 - user determined routes 


Random number seed used for stochastic attrition 


Width of search sector for unit i 
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ISPD 


ITEM 


ITIME 


ETRIT 


IUSTAT (I) 


LOA(I,J) 


LOST(I,J) 


LOT (I,J) 
MVTDIR (I) 
N (I) 

NA (I) 
NBU 


NF (I) 


NLOSC (I,J) 


NOI (I) 
NRU 
NT (1) 


OFL (I) 


It 


Input variable to denote user's desired speed 
for attackers movement 


l - 9 mph 
2 - 12 mph 
3 - 15 mph 
4 - 18 mph 


Input variable denoting number of time steps 
allowed for defender's move 


Current time, in seconds, of battle 
Input variable denoting whether attrition 
will be stochastic or deterministic 

0 - stochastic 

1 - deterministic 
Current status of unit it 

0 - unit alive, not firing 

1 - unit alive and firing 

2 - unit killed 

3 - unit moving 
The number of the ο. attacker of unit 1 


Denotes whether line-of-sight exists between 
unit i and j or not 


The number of the rn target of unit i 
Movement direction of unit i 

Number of nodes inputed by user for route i 
Number of attackers of unit i 

Number of defense units 


Number of time units i is allowed to fire at 
the same location 


Number of continuous time steps that line-of-sight 
does not exist between unit i and unit j 


Number of intervals in the ¡th route 
Number of attack units 
Number of targets of unit 1 


Force level of unit i during previous time step 
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E(I,J) 


PHH (I,J) 


PHM(I,J) 


PKH (I,J) 


PM 


POA(I,J) 


POL(I) 


PTT (I) 


RANGE 


Q(1,J) 


RF 


RMINTK 
'RMINTW 
RMXTK 
RMXTW 
ROT(1,J) 
SIZETK 
SIZETW 


TA (K) 


TF1(K) 


Probability of Js round hit by unit i in range 
band j 


Probability of a hit following a hit by unit 1 
in range band j 


Probability of a hit following a miss by unit i 
in range band j 


Probability of a kill given a hit by unit 1 
in range band j 


The proportion of time a moving unit is searching 
for targets 


The proportion of the a attacker of unit i 
allocated to fire on unit 1 


Percent of unit i lost during the current 
time step 


Proportion of surviving firepower allocated to 
the ith target if there are j targets available 


Current minimum distance between attackers and 
defenders 


Probability that unit j is not detected by unit i 
at current time 


Detection rate reduction factor for a firing unit 
(in comparison with non-firing unit) 


Minimum 


Minimum effective 


Maximum 


Maximum 


effective 


effective 


effective range for defending weapon 


range for attacking weapon 


range for defending weapon 


range for 


defending weapon 


system 
System 
System 


System 


th 


The range of the j target of unit i 


Size of attacking vehicle 


Size of defending vehicle 


Time to acquire a target for ο. 


type (k = 1,2) 


weapon system 


Time of flight to 1000m for gen 


type (k = 1,2) 


weapon system 
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—— [s  — ——4 


TF2 (K) 


ТЕЗ (К) 


ТН(К) 


TI (K) 


ТМ (К) 


TNKFR 
TOWFR 


TPOL(I) 


VISFR(I,J) 
VISFRA 
VISFRB 
E(I),Y(I) 
XA(I),YA(I) 
KIC(I,J) 
ETC(I,J) 
EL, YL 


KLOC (I,J) 
YLOC (I,J) 


Time of flight to 2000m for gh weapon system 
type (k = 1,2) 


Time of flight to 3000m for КС? 
type (k = 1,2) 


weapon system 
Time to fire a round following a hit for 
weapon system type k (К = 1,2) 


Time to fire first round after target has 
been acquired for weapon system type k (k = 1,2) 


Time to fire a round following a miss for 
weapon system type k (k = 1,2) 


Firing rate for attacking weapon system 
Firing rate for defending weapon system 


Total percentage of lost since battle began 
for unit i 


The fraction of unit i seen by unit j 
Fraction of unit A as seen by unit B 
Fraction of unit B as seen by unit A 
Coordinates of unit 1 


Coordinates of alternate position for 
defender i 


Coordinates of the jth interval endpoint of 


the route for unit i 


Distance added to previous interval endpoint 
for vehicle to move DST during a time step 


Coordinates of the ¿En node inputed by the 
user for the route of unit i 
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APPENDIX E 
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